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ABSTRACT 



To highlight in a qualitative manner the kinetics of an 
excimer laser, a simple computer model for calculating the 
species populations in a KrF laser cavity is developed; sub- 
sequently a computer program originally developed at the 
Naval Research Laboratory (NRL) is modified to calculate the 
population of the different electronic configurations of 
excited neon present in a XeF laser. When modified, the NRL 
program accounts for 135 reactions and requires 9.5 minutes 
of CPU time using the IBM 360/67. The populations obtained 
are applied to the calculation of the index of refraction 
in the laser cavity. The phase shift is determined per unit 
length for the two laser wavelengths; one laser wavelength is 
non-resonance with neon and the other is at resonance. Neon 
is the dominant specie relative to concentration and within 
the population distribution the neon ground state dominates 
by a factor of a million. The calculations show that the 
ground state neon determines the index of refraction; an 
exception occurs if the laser wavelength is near resonance 
to any of the transitions in the 3s -> 4p array. As long as 
the laser is operated away from the 3s -*■ 4p resonant wave- 
lengths, the phase shift will be negligibly small resulting 
in satisfactory beam quality; the preceding statement is 
valid only for the influence of neon. 
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I. INTRODUCTION 



The rare-gas monohalides, such as XeF, are simple diatomic 
molecules whose properties and emission spectra were essen- 
tially unknown as recently as five years ago. Now they are 
the active media for gas lasers that could provide overall 
electrical efficiency as high as 10 per cent in various wave- 
lengths , as depicted in Table I. 



Table I. Rare-gas Halide Laser Emission Wavelengths 
Peak wavelengths of most intense band (nm) 





Ne 


Ar 


Kr 


Xe 


F 


107 a 


193 


243 


351,353 


Cl 


b 


170 


222 


308 


Br 


b 


166 a ’ b 


206 a 


282 


I 


b 


b 


185 a ’ b 


252 a 



£ 

Has not demonstrated laser oscillation. 
b Predissociates , hence emission is weak or unobservable. 



The overall efficiency as explained by J . J. Ewing and C. A. 

Brau [1] is 

. . .a product of a quantum efficiency, an upper 
laser level production efficiency, an extraction 
efficiency, an efficiency for producing the initial 
excited states by the pump, and an energy coupling 
efficiency which describes the efficiency with which 
energy gets from a wall socket into the gas medium. 

The extraction efficiency is comprised of both 
spectroscopic and kinetic extraction efficiencies 
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depending on the ratio of net gain to loss and the 
ratio of rates of stimulated emission of upper laser 
levels to the quenching and spontaneous loss of 
upper levels. Certain rare gas halide lasers have 
potential for high power and high efficiency because 
they have high efficiency in each of the above 
mentioned elements. Depending on the application 
requirements, one may accept a lower efficiency of 
one of these elements if the compromise allows for 
some simplification of the laser. Typically this 
trade off usually involves coupling and excited state 
production efficiencies . 



Rare-gas monohalides belong to a larger class of mole- 
cules broadly known as excimers . An excimer is an atomic 
or molecular aggregate, composed of two atoms or molecules, 
that is bound in its excited state but is unstable or 
slightly bound in the ground state. 

The term "excimer" is standard chemical nomenclature 
for dimers which are bound in the excited state and free in 
the lower state. In standard nomenclature dimer means two 
of the same thing. Thus, Xe2*, A^* , etc. are excimers. 
Species that are heteronuclear such as LiXe* or KrF* also 
can have bound excited states and dissociative lower states. 
The accepted chemical nomenclature for such species is 
"exciplex" short for excited complex. Unfortunately, the 
laser community has not uniformly utilized good scientific 
English and typically calls all species with dissociative 
lower states, and even some like XeF with slightly bound 
lower states, "excimers". 

The excited state of the excimer species can radiate 
in a broad band that is red shifted from the wavelength of 
the parent atomic excitation. The size of the bandshift 



10 



depends on Che depth of the excited state well (typically 
l-5eV) and on the repulsion of the lower state at the energy 
minimum of the excimer. In Figure 1 are schematic potential 
energy curves for an excimer. 




Figure 1. Potential Energy of Excimers (Reproduced from 
Ewing [ 2 ] . 

Excimer emission spectra are useful for lasers because 

the lower level rapidly dissociates . Since the dissociation 
-12 

time, about 10 second, is much less than upper level 

-9 -6 

radiative lifetimes, 10 to 10 second, population inver- 
sions can be readily produced and maintained. 

The excimer laser or exciplex laser offers features 
which are essential for military applications [2-4]. These 
features are as follows: first, the potential for high 

efficiency; second, scalability to large volume and high 
power; and, third, transmission through the atmosphere 
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using either directly-generated photons or frequency-con- 
verted photons [5] . As a result of the attractive character- 
istics, considerable effort has been and is being directed 
toward identification of key technical issues and establish- 
ment of feasible solutions to the problems identified [6-11]. 
In addition to research directed toward the laser device 
itself, work is also being funded on related and relevant 
aspects, including e-beam technology, UV optics, and time 
dependent laser aerodynamics. For example, see Bennett, 
et al . [12] for work on UV optics. 

Systems studies have been conducted by Rocketdyne. Among 
the various aspects of large excimer lasers considered in the 
Rocketdyne reports was beam quality as influenced by non- 
uniformity of index of refraction. The adverse impact of 
poor beam quality upon the overall system was stressed. 

For military applications the optimum pulse duration 
should be a few microseconds [13-15]. The pulse duration 
for some military applications is in contrast to the desired 
pulse duration of a few nacoseconds for inertial confinement 
to achieve thermonuclear fusion [1]. Even if a pulse dura- 
tion of 10 microseconds is used, the motion of the waves, 
e.g., cathode shock wave and anode wave, in the gas will be 
very small, i.e., of the order of 

l - a x = (300 m/sec) (10 ^ sec) = 3 mm (1) 

where a is the acoustic wave speed and x is the pulse dura- 
tion. However, this fact does not suggest that a pulsed 
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excimer laser will be free of inhomogeneity problems; exten- 
sive experience with pulsed electrical lasers operating in 
the IR indicates that beam quality may be a problem. To 
achieve high average power, a large pulse repetition fre- 
quency, PRF , must be used which requires a flowing laser. 

The higher the PRF, the faster the gas flow is. 

Gas flow in a high power, high PRF, excimer laser will 
probably be subsonic. In subsonic flow [16—13] , disturbances 
can propagate upstream. Even if the upstream base flow is 
perfect, the gas moving out of the cavity after a pulse can 
cause An/n within the cavity; n is index of refraction, and 
An is spatial variation of index of refraction. A large 
flush factor eases the problem of a nonuniform base flow; 
however, high efficiency demands a low flush factor. 

Nonuniformity of the lasing medium may have adverse 
effects distinct from the deterioration of beam quality. 

The work of Nerheim and Chen [19] at JPL indicates different 
levels of performance depending on homogeneity of the gas . 
Pulses into a homogeneous medium resulted in one level of 
power output; pulses into an inhomogeneous medium yielded 
another level (lower naturally) of performance. Details of 
the physics involved in this observation were not unraveled. 
The work by Nerheim and Chen was with an atomic copper laser. 
Consequently the results are not directly transferable to 
rare gas -halide lasers. 

Despite the many unsolved problems regarding excimer 
lasers, its use as a potential weapon constitutes a 
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revolution greater than the leap from naval guns to guided 
missiles. Some of the conceivable missions for the excimer 
laser are ballistic missile defense (both ground-based silo 
defense and space-based destruction of ICBMs in the boost 
phase), space-based antisatellite defense, and antiship 
missile defense (ASMD) [20]. The latter has the most pro- 
mise of relatively near-term development. Also, the radia- 
tion from certain excimer lasers has good propagation in 
sea water. 

Briefly, the laser offers the advantage of depositing 
an extremely high energy on a small area of a target. This 
destructive energy is delivered at the speed of light (300 
million meters per second which is ideal for a defensive 
system against missiles); whereas, a conventional missile 
travels to its target at a speed of a few times the speed 
of sound (roughly 1500 meters per second) . 

As alluded to earlier, beam quality has an important 
impact upon the overall laser system. Improvements in beam 
quality can be used to decrease overall laser size, shorten 
time for melting a metal, or decrease laser run time to 
perform a certain task. A well know fact is that beam 
quality depends on favorable illumination of the exit 
aperture; favorable illumination implies near constant 
phase and, to a lesser degree, near constant amplitude. 
Variation of the index of refraction within the cavity 
causes spatial variation in phase at the laser exit. This 
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nonuniform phase causes the peak axis intensity, Iff. to 
decrease according to 

= exp [-(2T;<j) rms /A) 2 ] (2) 

o 

where X is wavelength. The rms phase, <J> , is obtained by 

squaring and averaging the phase shifts across the aperture. 
Equation (2) is the Strehle Ratio [21] . The spatial phase 
relation is 



p = / nds 



( 3 ) 



where ds is an incremental distance along a ray in the 
laser cavity. The initial efforts to limit phase distortion 
considered the index of refraction in terms of the Gladstone- 
Dale formulation 



. n = 1 + Gp (4) 

where p is gas density, and G is the Glads tone-Dale constant; 
hence there is concern for flow and acoustic effects. Hogge 
and Crow [22] state that flow uniformity and reduced acoustic 
effects have achieved Ap/p = 5x10 ^ , which is the required 
flow uniformity. 

It has been suggested, however, that the requirement for 
flow uniformity should be stated in terms of An/n due to 
nonuniform pumping in the cavity [23]. In a typical flow 
mixture, the constituents of the medium (NecXeiNF^) are 
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present in a fixed ratio (750:3:1) [24]. Even in a gas with 
fixed chemical composition, pressure and translational tem- 
perature, and hence, fixed Ap/p, An/n can vary according to 
the ionic and electronic excitation of the medium caused by 
pump ing . 

Consider the Ladenburg formulation of the Kramers -Kronig 
dispersion relation [25] 



n-1 = 



e 2 N 



N, 



27rm c 
e 



l j k 



Jsi 



A - A 



kj 



kj 



[1 - 



N kSj. 

Vk 



(5) 



e 

c 

m 



e 



kj 



kj 



N. 

J 



electron charge 
speed of light 
mass of electron 
total number 

wavelength for transition between k and j 

wavelength of interest at which n is derived 

is oscillator strength for transition from k to j 

th 

is degeneracy of j energy level 

th 

population in j energy level 



There are i distinct species indicating specific electronic 
or ionic excitation present in the cavity. Index j sums 
over lower levels and k over upper levels of possible con- 
figurations. Identify the mole fraction as 
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( 6 ) 



X 



l 




the fraction of the total molecular population represented 
by species i. Finally, define a quantity 



T 

Z 



2 inn c 
e 



E E 

j k 



2 2 
Xf . A 
Kl 



n 



2 2 
X ~ X kj 



f kj ^ 



k s 



P-1 



N i Sk 



(7) 



In terms of the quantities defined above 



n-l - N £ X 4 T t (S) 

£ 

or the index of refraction in the cavity is the sum of the 
contributions of all the species present in the cavity. 
Since N £ = NX £ , 

n-l = N 1 T 1 + N 2 T 2 + . . . (9) 

The partial density is given by = N^m^ where m ? is the 
molecular weight of species. In terms of partial 

density, equation (9) becomes 

n-l = p^T^/m^ + p 2 T 0 /m 2 + . . . (10) 

Introduce the mass fraction for the gas mixture = p /p, 
where p is the density of the gas mixture. The result is 

n-l = p[Y 1 T 1 /m 1 + Y^/im, + ...] (11) 
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Equation (11) shows that n-1 is proportional to p; however, 
equation (11) further demonstrates that n-1 depends on mass 
fraction and index of refraction of each species. Even with 
uniform flow, the spatial phase relation across the cavity 
will not be uniform. Accurate modeling of the phase effects 
requires knowledge of the species populations and the con- 
tribution to the index of refraction of the species in the 
cavity. This thesis will explore a model for determining 
the species populations in a laser cavity and apply a model 
originally developed at the Naval Research Laboratory (NRL) 
to calculate the populations of the different configurations 
of excited Neon present in a XeF laser. 

First, in section II, a review of electronic transition 
gas laser kinetics is discussed. Second, in section III, a 
simple model for an excimer (KrF) laser is developed. Third, 
in section IV, the computer program developed and written at 
NRL by Johnson, Palumbo, and Hunter [26] is discussed. 

Fourth, in section V, a formulation of input reactions to 
the NRL program to calculate the populations of excited 
neon configurations is developed. Fifth, in section VI, 
calculation of the index of refraction is performed. Sixth, 
in section VII, the calculation of phase shift is performed. 
Lastly, in section VIII, the results are discussed. 
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II. KINETICS 



Electronic transition gas lasers (i.e. exciplex or 
excimer) generally involve complex, branching chains of 
energy flow from initial pump excitation of background gas 
to eventual extraction of laser radiation. Figure 2 illus- 
trates major energy flow pathways in e-beam pumped 
Ne/Xe/NF^ mixtures. These energy chains include processes 
such as: energy deposition by electron ionization and 

excitation; energy transfer by various collisional 
mechanisms, including neutralization, displacement and 
"harpooning", leading to formation of the upper laser 
level; and collisional mechanisms which intercept excitation 
energy before it can reach the upper laser level, or which 
quench that level before stimulated emission occurs. Also, 
one must consider extraction of stimulated emission includ- 
ing absorptions by various processes in the gas and reso- 
nant reabsorption by the lower laser level (if there is 
one) , laser cavity oscillation, and coupling through output 
mirrors. A time-dependent model of the kinetics of such a 
system must define and update species densities for a large 
number--typically between 20 and 60--gas components, includ- 
ing ground states, several excited states and ions of atoms 
and molecules, plus electrons and photons. 

The time evolution of these components is followed by 
the solution of a set of coupled first order ordinary 
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Figure 2 
Note : 




Major energy flow pathways in e-beam pumped 
Ne/Xe/NF 3 mixtures (Reproduced from Huestis [27]). 

Energy input from the electrons comes in along 
the pump arrows. The species in each box, plus 
the component along the reaction arrow, yield 
the result in the next box. The wavelengths of 
the emitting species are indicated. They may be 
subject to quenching by NF 3 or Xe (not shown). 

NeF* is presumed to predissociate. 



20 



differential equations; one equation is required for each 
gas component. To the set of equations defining species 
densities are added further equations describing the evolu- 
tion of the photon field, the gas temperature, and (in some 
models) coupled electrical circuit parameters. The species 
equations can be generally described as equating the rate 
of change of species density to a sum of contributing for- 
mation terms less a sum of depletion terms as 

dN i V- V' 

air* E F ji - £ D i k < 12 > 

3 k 

Here is a particular species density, the F ^ ^ are the 

formation terms, and the the depletion terms. For 

example, consider a simple chemical reaction of the form 

A + S -+ C + D whose forward rate constant is k; if the 

species densities are expressed in cm , the units of k are 
3 

cm /sec. For this reaction, the rate equations defining 
species C and D would contain the formation term 

F - k[A] [B ] ; (13) 

at the same time, the equations defining species A and B 
would contain the depletion term 

D = k[A] [B] . (14) 

Recall that square brackets for a chemical symbols means 

3 

number density of the chemical with dimensions 1/cm . 
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Consider equation (12) again. Species i is formed in 
j different reactions. Likewise, species i is depleted in 
k different reactions. 
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III. A SIMPLE KINETICS MODEL FOR KrF 



In a thorough model of an excimer laser a large number 
of kinetic reactions are necessary; see, for instance, the 
185 reactions listed in the output of the NRL program for 
XeF (Appendix A) . The reactions are used to obtain the set 
of coupled first order differential equations. The solution 
to this set of equations is very involved and tedious 
because of the difficulty in handling the tangle of reactions. 
Therefore, only a few of the kinetic reactions in a KrF 
laser are considered. The object of the simple computer 
model was to obtain a qualitative insight to excimer laser 
kinetics . 

The model uses a He/Kr/F^ mixture (98.1/1.77/0.13) with 
He as diluent. Figure 3 illustrates a simplified e-beam 
pumped KrF laser geometry. 






Point of StlaulAted 
fiaiasion 



LAS£H CAVITY 




Figure 3. Example of e-beam pumped KrF laser geometry. 
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The first step in the kinetics chain leading to extrac- 
tion of laser radiation is the deposition of energy in the 
diluent gas by either photons or electrons. In this model 
electrons are used. The electrons ' ionize the Kr and 
dissociate the F 2 

k l + 

e' + Kr + l Kr + e + e' (15) 



where e' is a primary electron, e is a secondary electron 
resulting from ionization of Kr, and is the forward 
rate constant. The secondary electron has sufficient 
energy to dissociate and electronegatively attach to 
atomic F. 

k 2 

e + F 2 - F + F" (16) 

In addition dissociation of F 0 by the primary electrons 
occurs . 



e’ + F 2 - F + F~ (17) 

Having produced the ion pair states, recombination of 
positive rare gas ions and readily formed halide ions yields 
the excited state (KrF*) . 

+ - k4 

Kr + F + M - KrF* + M (18) 



The KrF* can produce a photon by stimulated emission. 



KrF* + hv 



KrF + 2hv 



(19) 
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Also, recombination of F can occur 



F + F+ M -*■ f 2 +m (20) 

Fluorescence decreases KrF* population by 

x 

KrF* - KrF + hv (21) 

where x is the radiative lifetime. 

Reactions (15) to (21) are the only reactions considered 
in this model. None of the quenching processes of KrF* 
have been considered, such as 

k 

KrF* + 2Kr Kr 2 F* + Kr (22) 

nor have many other reactions [2] and [28]. 

The reactions considered in this model are used to find 
the formation terms equation (13) and depletion terms 
equation (14). For example, examining reaction (15) the 
formation term for the concentration of e and Kr + is 
determined . 



Formation e = k^[e'][Kr] 


(23) 


Formation Kr + = k^[e'][Kr] 


(24) 



where [ ] denote concentration in specie number/cm . 

The depletion term for concentration of e from reaction (16) 
is 
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Depletion e = k 2 [e][F 2 ] 



(25) 



and for Kr + from reaction (18) 

Depletion Kr + = k^ [Kr + ] [ F~ ] [M] (26) 

Applying equation (12), that is, equating the rate of change 
of e concentration to the sum of contributing formation 
terms equation (23) less the sum of depletion terms equation 
(25) yields 



d[e] 

dt 



= k 1 [e*][Kr] - k 2 [e][F 2 l 



Similarly, combining equations (24) and (26) 

J - = '^[e’HKr] - k 4 [Kr + ][F‘][M] 



Continuing the process 



(27) 



(28) 



d[Fl 

dt 



d[F] 

dt 



= k 2 [e][F 2 ] + k 3 [e , ][F 2 ] - k g [F][M] (29) 

= k 2 [e][F 2 ] + k 3 [e'][F 2 ] - k^[Kr + ] [F~ ] [M] 

(30) 



djKrF*] = k ^ [Kr + ][F - ][M] _ k 5 [hv][KrF*] - [ Kr F* ] / x 

(31) 

^[ KrF] = k 5 [KrF*][hv] + [KrF*]/x (32) 



are obtained. 

In reaction (15) the e' is the e-beam input and there- 
fore is a function of time. The model assumes the 
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concentration of e' is initially 1.0E13 electrons/cm then 
becomes zero after about 177 nanoseconds. The secondary 
electron mean energy is 10 eV . 

The rate constants (k) for the reactions were taken 
from scientific articles [2, 29, and 30] or derived as in 
the case of . The values of the rate constants are listed 
in Table II. The radiative lifetime, t, is equal to 15 
nanoseconds from reference [ 2 ]. 



Table II. Rate Constants 






Rate Constant 


Value 


Reference 


k l 


5 . 0E-09 


3 

cm /sec 


29 


k 2 


2.0E-09 


3 

cm /sec 


30 


k 3 


5 . 0E-10 


3 

cm /sec 


30 


k 4 


1.0E-25 


cm /sec 


0 


k 5 


1.5E-05 


3 

cm /sec 


Derived in 
Appendix 3 


k 6 


1.0E-31 


6 , 

cm /sec 


2 


Equations (27) 


to (32) can 


be solved to 


calculate the 



concentrations as a function of time. Multiplying these 
equations by dt yields a form for numerical solution on the 
HP9830. Equations (27) to (32) become 



d [ e ] = (kj.fe'HKr] - k 2 [e][F 2 ])dt (33) 

d[Kr] = (k 1 [e'][Kr] - k^[Kr + ] [F~] [M] )dt (34) 
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d[F] - (k 2 [e][F 2 ] +k 3 [e'][F 2 ] -k 6 [F][M])dt (35) 

d[F“] = (k 2 [e][F 2 ] +k 3 [e'][F 2 ] - k 4 [Kr + ] [F~ ] [M] ) dt 

(36) 

d [ KrF* ] = (k 4 [Kr + ] [F~ ] [M] - k $ [hv] [KrF*] - [KrF*]/x)dt 

(37) 

d [KrF] = (kg [KrF* ] [hv] + [KrF*]/t)dt (38) 

respectively . 

The initial conditions for the set of equations are: 

[e] = 0 (39) 

[F] = 0 (40) 

[F~] = 0 (41) 

[Kr + ] = 0 (42) 

[KrF] = 0 (43) 

[KrF*] = 0 (44) 

Several constants are specified as follows: 

[hv] = 4.0E11 (calculated in Appendix C) (45) 

[F 2 ] = 4 . 0E16 'j (46) 

/ Assumed values for 

[Kr] = 5.4E17\ the He/Kr/F 2 mixture (47) 

( (93.1/1.77/0.13) at 

[He] = 3.0E19 1 STP . (48) 

Initial T = 0 

Final T = 250 nanoseconds 

Delta T = (Final T - Initial T)/256 

A computer program was written using equations (33) to 
(38) and the initial conditions and appropriate constants 
given by equations (39) to (43). Figures 4 and 5, which are 
plots obtained from the HP9330 program, are the concentrations 
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in specie number/cm versus time from 0 to 250 nanoseconds. 

The plots in Figures 4 and 5(b) indicate that some species 
densities rise extremely steeply in the first tenth of the 
simulation time; later in the pulse (after about 50 nano- 
seconds) species densities tend to flatten out. See Figure 4. 
Once the e-beam terminates, the decay time is approximately 
50 nanoseconds. In Figure 5(a) and (b) other species 
densities increase gradually and flatten out. The results 
of this model are in qualitative agreement with the NRL 
simulation contained in reference [26]. The NRL simulation 
limited the change in concentration to 5 percent per time- 
step to obtain sufficient accuracy [26]. The question of 
timesteps for the simple model was avoided; sufficiently 
small steps were used throughout the program. 

Table III lists the variables, symbols, and equations 
used in HP9830 computer program. 

Table IV lists the KrF Computer Program and Table V 
lists the output. 
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3.50E1 1 * 




(a) Concentration of secondary electrons 




(b) Concentration of Kr + 




(c) Concentration of KrF* 



Figure 4. 



Concentration of secondary electrons (e) , Kr + and 
KrF* , molecules /cm^ , as a function of time. 
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5.0E15 




(a) Concentration of F 




(b) Concentration of F 




Figure 5. 



(c) Concentration of KrF 

Concentration of F, F , and KrF 
as a function of time. 



molecules/cm' 
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Table III. Variables, Symbols, and Equations Used in 
HP9830 Computer Program for KrF Laser. 



Variables and Symbols Used in the Program 



HP9830 


Text 


Meaning or Value (for constants) 


Units 


El 


e* 


primary electron 


- 


E 


e 


secondary electron 


- 


F 


F 


Fluorine atom 


- 


G 


F 2 


Fluorine molecule 


- 


H 


F" 


Fluorine atom with extra electron 


- 


K 


Kr 


Krypton 


- 


L 


Kr + 


Krypton minus an electron 


- 


M 


M 


Third body 


- 


X 


KrF 


, Krypton Fluoride 


- 


Y 


KrF* 


Krypton Fluoride excited 


- 


P 


hv 


photon 


j oules 


K 1 


k l 


Rate constant 


3/ 

cm /sec 


K2 


k 2 


Rate constant 


cm^/sec 


K3 


k 3 


Rate constant 


cm^/sec 


K4 


k 4 


Rate constant 


6 i 

cm /sec 


K5 


k 5 


Rate constant 


3 

cm /sec 


K 6 


k 6 


Rate constant 


6 / 

cm /sec 


T1 


X 


Radiative lifetime 


sec 


I 


- 


Index; 1=1 for initial conditions 


- 


T 


t 


Time increment 


sec 
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Table III (Continued). Variables, Symbols, and Equations 
Used in HP9330 Computer Program for XrF Laser. 



Equations Used in the Program 



Line Number Equation Number 



860 


33 


370 


34 


880 


35 


890 


36 


900 


45 


910 


37 


920 


38 






X 
L U 
CO ~ 

x 

X * 



CO :t: 






: ! 




x 






f— i 




L3 ■¥ 










5— :t: 






= 




“7" 






il 




X + 
















*1— 1 * 




UJ * 










_J 










UJ -4: 


















s 


UJ ‘-4: 






j — 1 




LU 






— * 


Ui 


X 


: 




i 


bO 


X 






_l 


o 


h- * 


0— 




z 


UI 








II 


X 


X -4: 










: : 


UJ 








X -4 : 


_? 




-^-4 


u 


: * 


X 




X f— 




CO CO :4: 


= X 




: : 




LO X 


X 




- Vi 


o 


CM O ; t: 


X X = 




? — : ■ — : 


CO 


i 5 — 5 


x CD 




— • — * LU 


oo 


X 1 — # 


5—* 4: X 




; • l 




x 


X -• 




X W =—5 


X 


— 3 :t: 


lli i— m 




z z r- 


X, 


CD CD Q: 


X XX 




If II X 




UO UO LU -4: 


= -x - 




X 


0) 


CM CM : 4 : 


- l—« M- 




-! ! 1—i i i 




-—i O - 


IX X CO X 




X 


X 


j- X s— s -4 : X 


= 1 — X - • - 




X X 




-• - i— CO 


UJ i— * O CO Oj 




z = u i 


XI 


: *""! \jj M CO 


• : s : -iU 




r. 3-. > 


o 


CD 5 X X :4: CP. 


1 — x o 




? ? * : «! 1 




UO LiO 5— « O Z 


X X CO CM i-H 




t-J ■*““4 t- 5 r— 


WD 


CM CM X :4: X X 


Ll! = X H- 




: i • : • i X 


c 


• i £ i = W 


X LU LU X - ^ 




X X X s— « 


•H 


X X 5=4 :4: X X LU 


= = JPP O ^ Xi 




= = = =— i 


X 


- -• UJ O O h- 


CM CO j — i : — i i — * ~7~ X X 




II II II X 


cn 


r ” 7 ^ x :4: Ll. X O 


X» = X X X X o 




..... 


•H 


X* CD Cl. I— X 


CD CO *•-•* ; _= X 




-= — ! -= — * r, — S 


hJ 


in m x ; 4 ; co o x 


_J ^ X CO 




— * O J TT 




CM CM O *— • LU X 


CO X X j- X X 1 — x 




X X X f— X X 




i j t i O -4- 1 


X PLLl X X 






• 


xx x x x = 


jjJ X = ' « :_L. * f > r. 




- - - - II il 


> 


-• - X * X LLU 


X X X X X X X I — i — ? — z aa-. 




s 1 r-1 • ! ! |-;.j 1 ^- 


M 


i—l (— 1 , CU 2;- ,X X 


! 1 ! f O Q O i”i X* CO ! — 1 SI - 




r4 i-l H ^ V X 




CD CD X -4- X X :X ~ 


XXS— 1— XH-h-XX X X II 




• i i i i i • i = = 


<D 


iio i n o x x x 


i i i i X X X X X i — i x 1 1 — « X 




X X X -• -• 


c-H 


CM CM X -4' X X O O 


i » j— i — t J— j— X s X ! ? — ! | 


! D 




JO 


CO i_* L_s X X < O X 


X X X X X X X X h- f- >x 


\Tj 


II II II II X X 


cd 


-XXX -4- XXX 


XXXXr-XXXX 3 *— * X 


CM 


.•“. •. z z 


H 


X >•. JX CO X CO X 


1- X X *-i 


\ 


T—! I I * || IS 




i — i •' * * 5 O :-fc f— t 


x X »-* X 




^-4 CO 




X CD CD X X 


x x x x x >•: x x 5-h Gj .T-. 


CO 


X J- X X X X 




r- uo uo x -4* f— x : x x 


= = = = = = = = = = liO ! — CC' ! — CO 


i 


z z z z z z 




CM CM :::::: 


-T-, : : 


iT, 






! — • .■ • i 


1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — i — i — 1 — CO 


1— 1— 1— 1— 1— 1— 




X XXI — I — f — I — 1 — 1 — 


XXXXXXXXXX" X O X X« X x« 


■> — i 


X X X X X X 




X X X X X X X 


! i 1 i 1 ! f = ! * f * 1 } ! SI «ZD 1 * 1 CO X CO X 


II CO 


• — ! ! — i i — i . — s : — i i — ( 




X X X * ; — 1 * — } * — 1 ■ — * i — * 


X X X X X X X X X CM il X O * X *— • X 


— • II 


X X X X X X 




O ; ( ! — » X X X X X X 


x x x x x x x x x i — r* j x w s s-t p ^ 


s 

T 

p 

p 

p 

p 

p 

p 




XGPPXXXXXX 












<*H« iZ“j iZ“i i£;« ;3“t ! -Z~! i£l! ;Z“i iZ^i O Q iZ“i >£|i iZ“i i^ij iIZi 


CO X« X 1 X X' X '*u X ! 




X* X« X 5 X* X ! X X* X*' X 1 


•H* ♦ CM CO M" UO ! -£* p*- CO 'O' 1 X* ■*— * CM r- CO «7*‘ x> 


^ CM 


CO M" LiO *.jD P*- CO 




CM CO M“ UO CO P-- CO 1 J* 


i — t t— i -i — » i— ♦ -j~t -i — i t— » -r-H i— ( -r— t c*.j CM CM UO UO 1 UO UI f 


'.£• ‘O 


'•£* ’•£* '•£• *•£* ' £« X 
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X 



* 



4 



4 



4* 



: t- s ; 

0 8 i X 
4- 4 

■ : i ! 5 

•=— « !— < : f- 



4* 



4 4 s-i 
|_ Tf !_• 

TT” 

: : ! ■ : ! 



4 



i X *“* L2 X 



X 4 



U-: 

■ n 



I 


j— 






If 


X X 


: : 




s— 


X X 


■s—f 






S— 8 =— i 


i i 




X 


; — x 


X! 




X 


X 5 


a-. 




-•* 


1 — O 


: — ! 




f — i 


X X 


-T—i 




X 


CO 


!— i 




X z 


CO ■*- ■ 


X 




a*. 


X X 


>•. 




= CO 


i i 


; — s 




!i = 


: — i 


i i 


•:Ti 


CD si j — 


r— 


* 


CO 


X 1 — : 


X 


X 


T— ( 


: Si 


X 




-r— 4 


**• — 1 1 — 


■ — X 


CO 


CO 


UO X 






UO X 


X X 


X « 


CD 


Oi X 


- I — H- 


O S— 




X 


II >-« X 




a*. 


O I — 


UO X X 


o 


UO 


! 1-1 




o 


-r-J 


x — < 


: : : 


z z 




■’ — * UO ! I 






-r— J t— ( t—i \ [ 

l_i l_i i_i r - j x 



' — . — . — I — I — i — I — LU I — s— » 



!! OJ X X 
il If !! 



5 X X 

fi if 



X X X X 

XX X X 

CD X' UO X* 

jj-. Xi x ■>— < 
x r«- r- r-- 



x x X ! 
x x x : 

iH; -r-4 

OJ CO m 

r-- r- r- 



iX 

i i 



UO *-0 
CO CO 

r*- f-- 



X i i 

*-» «x 

x* x« 
uo co 
p- r- 



x* x x 



CD X 1 

CO x» 
r- co 



if II 

r- j x : 

x X* 1 
t— i C'-J 
CO CO ■ 



3 s 


4 4* 4 


If: 












■ ■ 


- — ! 5 . ! — ! . — ! 


UO -r— i 












X 


r : = j— i j i 


x \ — 












a-. 


' ! 4 : i ( L 1 


i X 












: ! 


X *“* X X 


. — 5 . — ! 












i — s 


4 : • 4 -r- 


3 i . 












• i 


o-i :_i oo tc 


3 i i i 












X 


X X X X 


J- 














I 4* 4* i 


:4: 4- 












■ — ! 


5- « C‘-i • !_! r_! 


? ? : ; 












J i 


: — » X « — 5 j — ; - 


. — ; = — t 












: i 


i i t- i i i i 


i 3 ■ i 












X 


X ^ L3 x 


X X 












a*. 


4 s-« 4 4= 


:f: 












3 X 


: ! . . i i i — ; 


. — r • — * 












: — i 


3 — * ijl : { -T— 1 


? i 1 












• — i 


! i : t : 1 1 ! l 


i 1 ' 1 












X 


3X *-* X X 


x 












a*. 


r hh 


4: 4: 




7T 








X 0 J 


t-h i ! cm '*“■* 


T UO 




:J* 








CO f — 


X X X X 


V 












CO 








~7~ 








1 x 


+ + + + 


4- 4* 




X 








t izr < 


1 ! I ! 1 ! 1 1 


1 ! i — 1 




X 








x - 


■ 1 i — 1 i 1 I — ( 


1 — 1 i — 1 




1 — 








X uo 


I I i i i i i *"“** 


1 i i i 












! — y—* 


X X X X — < 


j— X 




CD 






XI 




ii \\ if ii t_i 


ii ii 


S — CO 






7T 


X 1 


I 1 ! ! I ! I ! 


i ! I ! 


+ 


T— f 


* 


CD 


ij*. 


XXI — 


t-H -r*H -i— < t**H |{ 


r—t 


C'.j 






II 


II 


Z-- 1 — X 


+ + 4- 4* «- > 


4- 4- 


I— 


* 1 




i — i 


O «“* 


:■ — i 5 — i 


I — « l 1 j 1 3 ( !— 1 


• — * i — i 


if 




x 


t— t 


i — -=— • 


XXX 


i 1 * 1 1 i 1 i l ! 


i — i !—i 


0 i X X 


5_s 


O ! 1 


X X X X X X X 


0 — o-c 


1— 


i 1 


X 


X X X 


X* X UO CD «H i II 1 «X< X s 


CD X 1 


X« UO 


■II 1 


1 — 1 


0 J X5 


X U .5 uo 


‘.£i p - CO -‘J' 1 X 


?;*.j CO CO T3- Tj- Tj" UO 


CO CO CO 


CO CO CO CO X* 


IT. it, 


ij-. ij-r ij-i 


iT'« 


iT*« |J-, 
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Table IV (Continued) . Listing of the HP9830 KrF Program 



IT. Ill 









































iT. 




iT*. 




cj 




*j2* 








t-H 










































| — 




o 






*3 


Li”* 








t-H 


















r- 
































d* t-H 


-r-H 


OJ 




■ ! 






















t-H 
























! 




r— 




j*. \ 








: — i 




• 














*rH 


l -£* 


+ 


•T*. 


CO 




rr» 










i n 






**"! 




i 




rr. rr. 


rr. 


0 




i : 




~7~ 














t-H 


n^H 


LU 


-r-H 


-r-H 




d* 










d* 






LL 




CL 




s— 


n 


! 




*;.* 




o 














4* 


4* 


rj- 


+ 


4- 




I 


IT*. 


lT*« 


*21* 


do 


1 


t-H 


*2? 




IT’. 


~7~ 


*T-, 


3-. 3*. 








3*. 




i* ‘l 














LlI 


LU 


> 


LU LU 




LU 


d* 


d* 


t-H 


o i lu «:•:« r- 


5—4 


I — 


i— H 


o 


*21* *21* 


d* 


t-H 




*2*2* 




: 












•Hi 


*•3* 


Ti* 


LO 


CO 


t-H 




LiD 


1 


1 


1 


i 


Lf« 


1 


LiO 


z 




Z 








II 


• — 


H- 


1 — 1 






CD 


d* 


d* 


i^Z* 


II 


II 


II 


II 


’ll 


II 




• 


LU LU LU LU 


« 


LU 






1— 




1— LU CO 


»jj 


i — i 


:+: 






1 — I — 




il 


II 


ii 


{{ 


i — t 


r— * 


I — ! 


i — i 


1 — i 


* — i 


CD 


t-H 


LO 


01 it* 


t-H 


*i-H 


t-H 


i~« 


Ll- 


i 


r« 


~~i 


-J *— * 


■— « 




s — i 


1 — 


1 


~7~ CL 




i — i 


I — ! 


f — i 


I — ! 


-r-H 


-i — * 


t-H 


-i-H 


t-H 


t-H 


II 


II 


II 


II 


II 


II 


II 


I! 


! — 


CC* 


LL 


CO 


LL 


cc >c 




LL 


II 


o 




~7~ i — 1 !— 1 d* 


i— i 


T— I 


-r-H 


T— I 


-r-H 


1 — i 


! i 


1 1 


i l 


1 i 


K 2 


OJ 


-r-H 


t-H 


OJ 


CO 


rr 


U” 


«-£* 


o 


i — t 


~7~ 


t — i 


~7~ 


o »n 


»r 


o 


02* 


1 


LU LU LL LL 1— 


~7~ 


i i 




* 1 


i t 




LL 


f_2» 




2“ 


*2C 


j — 


! 








V 






i_n 


r—i 


5 1 


5—1 


; — i 


CO X 


j- 


Ll 


1 — 


;~i 




LL LL LL CO 


LU 


LL 


DI 


j 


- X - 
























































d' 


d* 


d* 


rri 


d* 




OJ 


CO 


in 


Li _i 

U” 


CD 




01* 






t-H 


u”. 1 ’ 


OJ CO H- IJ-I 


>3 


r£ *Ec* 


ij-. 


d* 


t-H 


«. *J 


CO 


^ IT* CO P- 


Li” 


CD 


r- 


CO 


CT*. 


Hi 


d* 


d* 


d* 






Hi 


d* 


iHi 


d* 


1^1 


t-H 


t-H 


t-H 


t-H 


t-H 


t-H 


t-H 


-r-H 


t-H t-H 


t-H 


OJ 


OJ 


OJ 


OJ 


O J O J O J O J CO 


iT. 


iT-. 


Cr* 


ij-. 


t-H 


**— • 


i-H 


-i-H 


t-H 


t-H 


-r-H 


-r-H 


t-H 


-r-H 


t-H 


t-H 


t-H 


t-H 


t-H 


t-H 


t-H 


t-H 


t-H 


t— H t-H 


t-H 


t-H 


t-H 


t-H 


t-H 


t-H t-H t-H -r-H 


t-H 
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Table IV (Continued) . Listing of the HP9830 KrF Program 



PROGRAM FOR COUPLED KINETIC EQUATIONS FOR THREE ELEMENT GAS LASER 



PT: 

U 






X 

-f- 

UJ 

;\,i 



! 

X 

3s 

3* 

3 













3 ; 


liO 












•30 


= 


t: 










3* 


rr— r 


: t : 












II 












=3 xT 




:i: 












X 












H !! 


X f\- 


* 












: 3 : 












-c— T T— r 


: — 1 












-w- X X 












X X T-4 


X x 


: t : 










i- 


: — i U "s UO X s 












X X X Aj CO *3 












=3 


1 1 \ 3 






j— 






3 ix X IX X UO 






•X 






iX* 


X 5 X 5 X 5 = 






X 






>3 


i — j .3 3* 3* 0 -1 


* 




x 






3« i — 3 3' 3* 






1 






, 


3 3 xi x! 11 






■3 






3* CO 


i — « 3' 3* 3- ; — 






X 








jZ^ > » 3 


: t : 




X 






H IS 


3 AJ — * X 






X 








X 3 








X 




■rS T-i 


H II It 3 






* 






•w- *w 


— * 0 J X 


* 










3 3 


h- X X X X 






1— 


UO 




-i- 




:t: 






x 




X 




:f- 




X 


a-. 




3* 








■x 






3s 






3* 


X CO 






3- 






0"i 


1 — i — X 


r. 




3 




* 


CO 


X 5— < o 


CO CO 


3* 






U-. 


X o o 


x 


~7~ 


• 






LU 


X X 




o 


3 ‘TT 






X X 
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Table V. Output of the 1IP9830 KrF Program 
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Table V (Continued) . Output of the HP9830 KrF Program 
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Table V (Continued) . Output of the HP9830 KrF Program 



IV. 



XeF LASER COMPUTER PROGRAKFROI1XRL 



In order to calculate the excited states of neon in a XeF 
laser a complex kinetics program is required. A kinetics 
scheme for the XeF laser developed at the Naval Research Lab 
was used. It was written by Dr. Louis J. Palumbo [26 and 31] 
who assisted greatly in modifying it for use on the IBM 360/67. 

The computer model consists of two parts: first, an e- 

beam deposition section, in which the kinetic energy of near 
relativistic electrons is used to ionize and excite the laser 
gas, and, second, a chemical kinetics section, in which the 
energy deposited in the gas cascades down by various processes 
to the upper laser level and is extracted as laser radiation. 

One input to the NRL computer model is a list of reactions, 
e.g. in the format of equations (16) to (21). Also an input 
is a list of reaction rate coefficients. The computer program 
sorts the reactions and forms formation terms, F, and 
depletion terras, D. Formation and depletion terms were 
discussed earlier in connection with equations (13) and (14) . 
Continuing, the computer program selects the appropriate F 
and D terms to form the rate equation of equation (12). Next, 
the computer program multiplies the rate equation by a variable 
times tep. The same approach was used in the simple computer 
program; see equations (33) to (38) . 
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An advantage of the NRL program is the ability to change 
the number of reactions. The computer program does not change; 
only the input changes. The differential equations are 
generated internal to the program. This method makes it much 
easier to alter the set of equations and greatly reduces errors 
attendant to such alterations which arise with ease in the 
complex tangle of reactions. The detailed explanation of the 
computer model is contained in references [26] and [31]. 
Appendix A gives a list of reactions used to generate the 
population distribution within neon energy levels. 
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V . CALCULATION OF NEON EXCITED 



STATE POPULATIONS 



The excited states of neon are represented in Figure 6. 
The electronic configurations of interest to this study are 
the 3s, 3p, 4s, 3d, 4p, 5s, 4d, 4f, 5p, 5d, and 6s. 

In order to calculate the populations of these particular 
configurations, one needs to know the reactions involved and 
the rate constants. For example, 

k i 

Ne + e-beam -*■ Ne(3s) 
k 2 

Ne + e-beam Ne(3p) 



Ne 4- e-beam 



’ll 



Ne (6s) 



These reactions could then be input to the NHL XeF pro- 
gram to calculate the populations of the excited configura- 
tions of neon. After, an extensive literature search (Optics 
and Spectroscopy, The Journal of Chemical Physics, Applied 
Physics Letters, Physical Review A, Journal of the Optical 
Society of America, Chemical Abstracts, IEEE Journal of 
Quantum Electronics, etc.) rate constants for the individual 
excited configurations i.e. 3s, 3p , ... were not found. 

Therefore, it was necessary to make a number of assump- 
tions. The procedure followed was to extract all reactions 
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Figure 6 




Energy-Level diagram of Neon (Reproduced from 
Pressley [ 32] . 



/. Q 



from the NRL program involving the neon excited state (Ne*) 
Then, these reactions were subdivided into two groups, the 
Ne(3s) and all higher configurations, which are denoted by 
NE*1 and NE*Z respectively. The population distribution 
within NE*Z was calculated assuming a Boltzmann distribution 
The list of reactions which were used as modified input 
to the NRL program are listed in Table VI in the following 
format : 



Ne 


+ 


e 1 


- Ne* 


+ 


e’ 


Rate Constant 
1.20E-13 


(49) 


Ne 


+ 


e 1 


- Ne*l 


+ 


e' 


0.925E-18 


(50) 


Ne 


+ 


e ’ 


- Ne*Z 


+ 


e’ 


0 . 275E-18 


(51) 



To compare equations (4-9) to (51) with the first three 
reactions of Table VI note that e' is EBEM in FORTRAN nota- 
tion. The other symbols are obvious. The first reaction 

yielding Ne* is divided into the Ne*l and Ne*Z reactions. 

2 

A computer run with the current density equal to 10 A/ cm 

and a computer run with the current density equal to 11 
2 

A/ cm calculated the populations for Ne*l (3s) and Ne*Z 
(all levels greater than 3s) . 

The rate constants for the reactions were determined by 
taking a fraction of the total rate constant. For example, 
in Table VI, for the second set of reactions, the rate 
constant for NE*1 is 0.3 and NE*Z is 0.2 of the rate con- 
stant for NE*. These fractions of the rate constant for 
NE* are noted in parentheses following the reactions. The 
choices of the fractional values were arbitrary. 
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TABLE VI. Input Reactions to the NRL Program. 



Reaction Rate Constant 



NE 


+ 


EBEM 






NE* 


+ 


EBEM 








1.20 E-18 


NE 


+ 


EBEM 




-+■ 


NE*1 


+ 


EBEM 








0 . 925E-18 


NE 


+ 


EBEM 




-V 


NE*Z 


+ 


E3EM 








0 . 275E-18 


NE* 


+ 


T? 




-► 


NE 


+ 


E- 








2.00 E-10 


NE*1 


+ 


E- 




-> 


NE 


+ 


E- 






(.3) 


1.60 E-10 


NE*Z 


+ 


T7_ 




-V 


NE 


+ 


E- 






(.2) 


0.40 E-10 


NE2+ 


+ 


E- 




-> 


NE* 


+ 


NE 








3.50 S-08 


NE2+ 


+ 


E- 




-+■ 


NE*1 


+ 


NE 






(.9) 


3.15 E-08 


NE2+ 


+ 


E- 






NE*Z 


+ 


NE 






(.1) 


0.35 S-08 


NE3+ 


+ 


E- 




-+ 


NE* 


+ 


NE 


+ 


NE 




3.50 E-08 


NE3+ 


+ 


J-l 




->■ 


NE*1 


+ 


NE 


+ 


NE 


(.95) 


3 . 325E-08 


NE3+ 


+ 


E- 






NE*Z 


+ 


NE 


+ 


NE ' 


(.05) 


0.175E-08 


XENE+ 


+ 


E- 




-> 


NE* 


+ 


XE 








6.00 E-08 


XENE+ 


+ 


E- 






NE* 1 


+ 


XE 






(.6) 


3.60 S-08 


XENE+ 


+ 


E- 






NE*Z 


+ 


XE 






(.4) 


2.40 E-08 


NE* 


+ 


NE + 


NE 




NE2* 


+ 


NE 








5.0 


E-34 


NE*1 


+ 


NS + 


NE 


-y 


NE2* 


+ 


NE 






(.8) 


4.0 


E-34 


NE*Z 


+ 


NE + 


NE 


-y 


NE2* 


+ 


NE 






(.2) 


1.0 


E-34 


NE* 


+ 


XE 




- 


XE+ 


+ 


E- 


+ 


NE 




7.40 E- 11 


NE*1 


+ 


XE 




-y 


XE+ 


+ 


E- 


+ 


NE 


(.6) 


4.44 E- 11 


NE*Z 


+ 


XE 




-y 


XE+ 


+ 


E- 


+ 


NE 


(.4) 


2.96 E- 11 


NE* 


+ 


XE 




-y 


NEXE+ 


+ 


E- 








2.20 E-ll 


NE*1 


+ 


XE 




-y 


NEXE+ 


+ 


E- 






(.6) 


1.32 E-ll 


NE*Z 


+ 


XE 




-y 


NEXE+ 


+ 


E- 






(.4) 


0.88 E-ll 
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TABLE VI (Continued) . 



Input Reactions to the NRL Program. 



* 



Reaction Rate Constant 



NE* 


+ 


F2 


-y 


NEF* 


+ 


F 








4.10 


E- 10 


NE*1 


+ 


F2 




NEF* 


+ 


F 






(.3) 


3.28 


E-10 


NE*Z 


+ 


F2 




NEF* 


+ 


F 






(.2) 


0.62 


E-10 


NE* 


+ 


NF3 




NEF* 


+ 


NF2 








1.05 


E-10 


NE*1 


+ 


NF3 


-V 


NEF* 


+ 


NF2 






(.8) 


0.84 


E-10 


NE*Z 


+ 


NF3 


-> 


NEF* 


+ 


NF2 






(.2) 


0.21 


E-10 


NE* 


+ 


NF2 


-h 


NEF* 


+ 


NF 








1.0 


E-10 


NE*1 


+ 


NF2 


-+■ 


NEF* 


+ 


NF 






(.3) 


0.8 


2-10 


NE*Z 


+ 


NF2 


-y 


NEF* 


+ 


NF 






(.2) 


0.2 


E-10 


NE* 


+ 


NF 


-y 


NEF* 


+ 


N 








1.0 


E-10 


NE*1 


+ 


NF 


-y 


NEF* 


+ 


N 






(.8) 


0.8 


E-10 


NE*Z 


+ 


NF 


-> 


NEF* 


+ 


N 






(.2) 


0.2 


E-10 


NE2* 


+ 


HN352L 


-y 


NE* 


+ 


NE 


+ 


ABSL 




1.0 


E-18P 


NE2* 


+ 


HN352L 


-y 


NE* 


+ 


NE 


+ 


ABSL 


(.9) 


0.9 


E-18P 


NE2* 


+ 


HN352L 


-y 


NE* 


4* 


NE 


+ 


ABSL 


(.1) 


0.1 


E-18P 



* 

Use of capital letters more closely represents computer 
printing. 
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However, the first reactions listed in Table VI do not 
have fractional values in parentheses, as the reaction rate 
constants were determined using reference [9]. The plot of 
"Excitation of Neon by Electrons" from page 467 of refer- 
ence [9] was extrapolated to 10 keV, and the value for the 
2p state was assumed to represent all excited configura- 
tions above the 3s (NE*1) . The difference between the rate 
constant for NE* and the 2p-^ value yielded the rate constant 
for NE*1. 
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VI . CALCULATION OF THE INDEX OF REFRACTION 



The values (n-l)/N where n is the index of refraction 
and N is the species population, have been calculated Refs. 
[33], [34] and [35]. Figure 7 lists the values of (n-l)/N. 
References [34] and [35] contain the details of calculating 
these values. 

The index of refraction (n-1) is equal to 



n- 1 = 



(Ijlzil) N + N(3s) + N(Z) 



(52) 



ground 

state 



Ne*l 



Ne*Z 



Using the populations calculated in the NRL model for N, 
N(3s), and N(Z) equation (52) can be solved. However, the 
populations of individual electronic configurations must be 
determined. 

One feature of excimer lasers, in contrast to the 
vibrational transition lasers of the infrared region, is 
the extent of randomness of populations. In the C0 9 gas 
dynamic laser, mode temperatures can be assigned to the 
various vibrational modes. Due to the short span of excita- 
tion and power extraction in an excimer laser, the various 
modes do not have an opportunity to equilibrate. This fact 
is the reason extensive computer programs are necessary for 
kinetics. An assumption of a Boltzmann distribution is 
unfounded. However, in view of the sparsity of reaction 
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Figure 7. Matrix of significant configuration arrays. The number in the blocks is 
(n-l)/N cm multiplied by 1C)27 (Reproduced from Fuhs, Etchechury, and 
Cole [33]). 7 




rate coefficients for excitation to specific configurations 
a Boltzmann distribution is assumed. Assuming a Boltzmann 
relation for the upper levels the populations of the 3p , 

4s, 3d, ... 6s can be calculated. 

N(Z) = N(3p) + N (4s) + N (3d) + ... + N(6s) (53) 



N(Z) = N(3p) [1 + 



g 



g 



-< E 4s- E 3p> 



4s 



kT 



3p 



exp 



‘3d 

: 3p 



exp 



- (E 3 d -E 3p ) 
kT 



..] ( 54 ) 



where g denotes the degeneracy of the configuration which is 
equal to the sum of the 2J + 1 terms. E is the energy sub- 
scripted by the configuration. 

Solving for N(3p) 



N(3p).= 



N(Z) 



1 + 



^4s 

s 3p 



exp 



-< E 4s- S 3p> 

kT 



+ 



(55) 



All the quantities on the right side of equation (55) 
are known or can be calculated. Table VII lists the values 
of the degeneracies and energies for the configurations. 
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TABLE VII. Degeneracy 
conf igurat 


and energy values 
ions . 


for the 


Electronic 

Configuration 


SC2 2J +1) 


E (cm -1 ) 


3s 


12 


134674.6 


3p 


36 


150197.8 


4s 


12 


158950.1 


3d 


73 


161874.4 


4p 


36 


163180.3 


5s 


12 


166123.2 


4d 


73 


166500 


4f 


84 


166562 


5p 


36 


167812.5 


5d 


73 


168000 . 


6s 


12 


168750 



The temperature T is calculated as follows using the 
Boltzmann relation 



H(3s) _ £3 s 
NWX g 2 p p 

Solving for T 



-(E- -E-, ) 
v 3s 3p 



•NT2pT 



’3s 



kT 



< E 3s -E 2p> 

kT 



( 56 ) 



T 



-(E 



3s 



JLiad. 



k Jin 



r N(3s) 

l NT2 pT 



S 



1^2 

3s' 



(57) 



Substituting the value of T into equation (55) yields a 
value for II(3p). To calculate N(4s) one uses 



N (4s) = 



N(3p) exp 



-^4s- E 3p> 

kT 



s 3p 



(58) 



Similarly the populations of the other configurations 
can be calculated i.e. IT ( 3d) , N(4p), etc. Once the popula- 
tions of the configurations have been found it is a simple 
calculation to obtain n-1 applying equation (52) as 



n-1 n-1 n-1 

Ue(3p) ’ He (4s) ’ He (5s) 

are given along the bottom of Figure 7. 

A computer program to calculate the index of refraction 
(n-1) using equations (52'J , (55), (57), and (58) was 

written. 

In Figure 7 the 3s -*■ 4p transition contains an R*. R* 
indicates that resonance may occur and can have a large 
effect on the index of refraction. The spikes in Figure 8 
are located at the resonant wavelengths. Resonance is 
discussed in reference [36]. 
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The plots in Figure 9 are of the calculated indices of 
refraction (n-1) versus time (nanoseconds). Figure 9(a) 
is for a laser wavelength away from 3s -*■ 4p (no resonance) 
for both low and high power. The current density for low 
power is 10 A/cm , and the high power is 11 A/ cm . For a 
nonresonant wavelength, the change in n-1 is in the fourth 
decimal place. 

The resonance effect was taken into account in a second 
computer run; see Figure 9(b). For a resonant laser wave- 
length, the contribution to the index of refraction due to 
Ne(3s) was calculated assuming (n- 1) /Tie (3s) equal to 
-5.0X08; the corresponding value for Xe(4p) was (n-l)/Ne(4p) 
equal to 1.667E08. Figure 9(b) is for a laser wavelength 
near the 3s -*■ 4p resonance for low and high power. The 
resonance has a marked effect on the index of refraction 
in comparison to Figure 9(a) as evidenced by the ordinates 
of the plots. Also, Figure 9 reveals that increasing the 
power decreases the index of refraction. 

The right hand ordinate scale for Figure 9(a) is the 
ratio (n-1) / (n-1) , where (n-l) Q is the value at zero time. 

The scale emphasizes the fact that the relative changes are 
in the fourth decimal. 

Table VIII lists the variables, symbols, and equations 
used in the HP9830 program to calculate n-1. Table IX lists 
the HP9330 program used to calculate n-1 and Tables X and XI 
list sample output for low and high power respectively. 
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(a) Index of refraction (n-1) at laser wavelength away from 3s 4p 
resonance for high and low power. 




(B) Index of refraction (n-1) at laser wavelength near 3s 4p resonance 
for high and low power. 



Figure 9. Calculated index of refraction (n-1) as a function 
of time. 
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(T-u)/(T-u) 



TABLE VIII. 


Variables, Symbols, and Equations Used in 
HP9830 Computer Program for Calculation of 
n-1. 




Variables and Symbols Used in the Program 



HP9830 Text Meaning or Value (for constants) Units 



G g 


degeneracy (I 2J+1) 


E E 


energy cm ^ 


Cl c 


speed of light cm/sec 


K2 k 


3oltzmann constatit eV/°K 



M n-l/Ne() Contribution of excited state 

K 2p,3s,.. Index of the levels 1 to 12 



K() T 


Temperature degrees 

Kelvin 


J 


Time increment index 


Line Number 


Equations Used in the Program 

Equation Number 


320 


57 


930 


56 


1030 


54 


1060 


55 


1130 


53 


1350 


52 
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TA3LL IX. Listing of the HP9330 Index of Refraction (n-1) 
Program. 



10 PRINT " LOW POWER 1 0P/CM t 2" 

20 PRINT 

100 DFiTEi 0? 9. 716E+19? IE-50? IE-50 

1 1 O DPT fl 5 01 ? 9 » 7 1 6 E + 1 9 ? 2 . 8 5 E + 1 3 ? 1 . lb b E + 1 3 
120 D P T P 1 01 0 ? 9 . 7 1 6 E + 1 9 ? 5 . 3 1 4 b. + 1 3 ? 2 . 8 0 2 E + 1 ’J 
130 DPTP 150? 9. 716E+19? 7. 709E+13? 4. 034E+13 
140 DPTP 2 0 0 .? 9 . 7 1 6 E + 1 9 ? 8 . 8 6 8 E + 1 3 ? 4 . 8 4 3 E + 1 3 
150 DPTP 250? 9. 71SE+19? 9. 572E+13? 5. 35E+13 
160 DPTP 3 00 ? 9 . 7 1 6 E + 1 9 ? 9 . 9 9 8 E + 1 3 ? 5 . b 6 1 E + 1 3 
170 DPTP 3 5 O ? 9 . 7 1 6 E + 1 9 ? 1 . 0 2 b t + 1 4 ? b . 3 5 1 E + 1 3 
180 DPTP 400? 9. 716E+19? 1 . 041E+14? 5. 966E+13 
190 DPTP 450? 9. 716E+19? 1 . 051E+14? 6. 035E+13 
200 DPT P 500 ? 9 . 7 1 6E+ 19? 1 . 057E+ 1 4 ? 6 . 077E+ 1 3 
2 10 DPTP 550 ? 9 . 7 1 6E + 19? 1 . 06E+ 1 4 ? 6 . 1 03E+ 1 3 
220 DPTP 600? 9. 716E+19? 1 . 0S2E+14? 6. 1 18E+13 
230 DPTP 6 5 0 ? 9 . 7 1 6 E + 1 9 ? 1 . 0 6 3 E + 1 4 ? 6 . 1 2 7 E + 1 3 
240 DPTP 700? 9. 716E+19? 1 . 064E+14? 6. 132E+13 
250 DPTP 750 ? 9 . 7 1 6E+ 1 9 ? 1 . 065E+ 1 4 ? 6 . 1 35E *• 1 3 
2 60 DPTP 3 O0? 9 . 7 1 b E + 1 9 ? 1 □ 0 6 5 E + 1 4 ? b . 1 3 7 E + 1 3 
270 DPTP 8 5 O ? 9 . 7 1 6 E + 1 9 ? 1.065b+14?S. 13 8 E + 1 3 
280 DPTP 9O0 ? 9 . 7 1 6E + 1 9 ? 1 . 066E+ 1 4 ? 6 . 1 39E+ 1 3 
290 DPTP 950? 9. 716E+19? 1 . 066E+14? 6. 14E+13 
3 0 0 DPTP 1 0 O O ? 9 . 7 1 6 E + 1 9 ? 1 . 0 b 6 E + 1 4 ? 6 „ 1 4 E + 1 3 
310 DPTP 1 0 5 0 ? 9 . 7 1 6 E + 1 9 ? 1 . 0 6 6 E + 1 4 ? 6 ,, 1 4 E + 1 3 
320 DPTP 1 1 O 0 ? 9 . 7 1 6 E + 1 9 ? 1 . 0 b 7 E + 1 4 ? 6 . 1 4 E + 1 3 
3 30 DPTP 1 1 5 0 ? 9 . 7 1 6 E + 1 9 ? 1 . 0 6 7 E + 1 4 ? 6 . 1 4 1 E + 1 3 
340 DPTP 1 2 0 0 ? 9 . 7 1 6 E + 1 9 ? 1 . 0 6 7 E + 1 4 ? 6 . 1 4 1 E + 1 3 
350 DPTP 1 2 5 0 ? 9 . 7 1 6 E + 1 9 ? 1.0S7E+14?b. 1 4 1 E + 1 3 
360 DPTP 1300? 9 . 7 1 6 E + 1 9 ? 1 . 0 6 8 E + 1 4 ? 6 . 1 4 2 E + 1 3 
370 DPTP 1 3 5 0 ? 9 . 7 1 6 E + 1 9 ? 1 . 0 6 3 E + 1 4 ? 6" . 1 4 2 E + 1 3 
330 DPTP 1 4O0 ? 9 . 7 1 6E+ 19? 1 . 069E+ 1 4 ? 6 . 1 43E+ 1 3 
390 DPTP 1 4 5 0 ? 9 . 7 1 6 E + 1 9 ? 1 . 01 6 9 E + 1 4 ? b . 1 4 4 E + 1 3 
400 DPTP 1 5 0i 0i ? 9 . 7 1 6 E + 1 9 ? 1 . 0 7 E + 1 4 ? 6 . 1 4 b E + 1 3 
410 DPTP 1 550 ? 9 . 7 1 6E+ 1 9 ? 1 . 07E+ 1 4 ? 6 . i 46E+ 1 3 
4 2 0 DPTP 1 6 0 0 i ? 9 . 7 1*6 E + 1 9 ? 1 „ 0 7 1 E + 1 4 ? 6' . 1 4 7 E + 1 3 
430 DPTP 1650? 9. 716E+19? 1 . 072E+14? 6. 148E+13 
4 4 01 DPTP 1 7 01 01 ? 9 . 7 1 6 E + 1 9 ? 1 ,, 0 7 3 E + 1 4 ? 6 . 1 4 9 !• + 1 3 
450 DPTP 1750? 9. 716E+19? 1 . 074E+14? 6. 15E+13 
4601 DPTP 1 800 ? 9 . 7 1 6E+ 19? 1 . 074E+ 1 4 ? 6 » 1 52E+ 1 3 
470 DPTP 1 850 ? 9 . 7 1 6E+ 19? 1 . 075E+ 1 4 ? 6 . 1 53E+ 1 3 
480 DPTP 1 900 ? 9 . 7 1 6E+ 19? 1 . 077E+ 1 4 ? 6. 1 55E+ 13 
4 9 01 DPTP 1 9 5 0 ? 9 7 1 6 E + 1 9 ? 1 . 0 7 3 E + 1 4 ? 6 . 1 5 7 E + 1 3 
500 DPTP 2000? 9. 716E+19? 1 . 079E+14? 6. 159E+13 
510 DIM TC 41 1 ? KC 4 1 ]?SC41 1 ? ZE 4 1 J ? ME 1 3 3 ? G C 13] 

515 DIM RE 4 1 ? 1 2 3 ? LE 4 1 ] ? HE 4 1 3 ? EE 1 3 ] ? PE 4 1 ? 1 2 11 ? NC 4 1 ? 1 2 ] 
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TABLE IX (Continued) . Listing of the I1P9830 Index of 
Refraction (n-1) Program. 



518 Y=0 " 

520 GC 1 3=1 

530 GC 2 3=GC 4 3 = GC 7 3=GC 12 3=12 
540 G C 3 3 = G [ 6 3 = G C 1 0 3 = 3 6 
550 GC 5 3=GC 8 3=’GC 11 3 = 73 

5 6 y G C 9 3 = « 4 
570 EC 1 3 = 0 

530 EC 2 3=134674. 61 
590 EC 3 3=150197. 31 
630 EC 4 3=153950. 1 
610 EC 5 3=161374. 4 
620 EC 6 3= 1 63 1 30 . 29 
630 EC 7 3=1 66 1 23 . 23 

6 4 0i E C 3 3 = 1 6 6 5 Gi 0 
650 EC 9 3= 166562 
6601 EC 1 01 3= 1 67 y 1 2 . 5 
670 EC 1 1 3 = 163000 

6 3 0i EC 12 3=16 y 7 5 6i 

690 K2=3. 6171E-05 

69 1 C 1 =2 . 997925E+ 1 0 

692 H=4. 1356E-15 

693 C2=C1*H 
700 C= 1 E-27 
705 Y=Y+ 1 

710 FOR J=1 TO 41 

720 READ TC J 3? NC J? 1 3 j NC J j 2 3 ? Z C J 3 

730 NEXT J 

740 FOR J=1 TO 41 

770 NEXT J 

790 FOR J=1 TO 41 

8 0 0 Q = N C J » 2 3 * G C 1 3 / ( N C J > 1 3 * G C 2 3 > 

310 Q=LOGQ 

8 2 0 K C J ]=-< (EC 2 3— EC 1 3)* C 2 ) / < Q * K 2 > 

350 NEXT J 

380 FOR J=1 TO 41 

900 FOR K=4 TO 12 

9 1 0 Q 1 = < ( E C K 3-EC 3 3) * C 2 > / ( K 2 * K C J 3 > 

920 Q1=EXP< -Q1 > 

930 AC J j K 3=GC K 3*Q 1 -"GC 3 3 

960 NEXT K 

980 NEXT J 

990 FOR J=1 TO 41 

1000 3=1 

1020 FOR K=4 TO 12 

1030 S=S+AC J ? K 3 

1040 NEXT K 

1050 SC J 3=3 

1060 NC Jj3 3=ZC J3/3C J3 

1090 NEXT J 
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TABLI IX (ConCinued) . Listing of the HP9830 Index of 
Refraction (n-1) Program. 



1100 


FOR J=1 TO 


41 


1120 


FOR K=4 TO 


12 


1130 


NC Jj K ]=NC J j 


3 1*RC J i> K 1 


1 160 


NEXT K 




1170 


NEXT J 




1210 


MC 1 1=558 




1220 


MC 2 ] = -87444 




1230 


MC 3 ] = -59674 




1240 


MC 4 ] = -58722 




1250 


MC 5 1=31324 




1 260 


MC 6 l=-6 1 758 




1270 


MC 7 l = -549 1 8 




1 280 


MC 81 = 27817 




1 290 


MC 9 1=0 




1 300 


MC 10 l=-3243 


5 


1 3 1 0 


MC 1 1 1=48799 




1320 


MC 121=84238 




1330 


FOR J=1 TO 


41 


1340 


R=0 




1350 


FOR K = 1 TO 


12 


1360 


RC J ? K 1=MC K 1 


*NC J j K 1*C 


1370 


R=R + RC J? K 1 




1400 


NEXT K 





1412 IF Y= 1 THEN 1418 
1414 HC J ]=R 
1416 GOTO 1420 
1-418 LC J ] = R 

1420 PRINT "N - 1 I S " R " R T T I ME " TC J ] " NSEC 
1440 NEXT J 
1450 PRINT 

1460 PRINT " HIGH POWER 1 1 R/ CMt2 

1465 PRINT 

1468 IF Y >= 2 THEN 1500 
1470 GOTO 705 

1 500 SC RLE 0 ’ 2000 j -2E-05 ? 6E-05 

1510 XRXIS -2E-05) 100 

1526 Y FlX IS 0 j IE-05 

1530 FOR J=1 TO 41 

1540 PLOT TC J ] j LC J ] 

1550 NEXT J 
1560 PEN 

1570 FOR J=1 TO 41 
1580 PLOT TCJ'JjHCJ] 

1590 NEXT J 
1600 PEN 
1610 STOP 
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TABL3 



IX (Continued) . Listing of the HP9830 Index of 
Refraction (n-1) Program. 



5080 


PRINT 


" DflTfl 


5 0 1 0 


DfiTfi 


0 ? 9 . 7 1 1 


5020 


DflTfl 


50? 9.71 


5030 


DflTH 


1 00 3 9 , 1 


5040 


DflTfl 


1 5 U ? 9 . i 


5050 


DflTfl 


200 ? 9. 1 


• j 0 b 0 


DflTfl 


250? 9. 1 


5070 


DflTfl 


300 ? 9 . 1 


5080 


DflTfl 


3 5 0 ? 9 . i 


5090 


DflTfl 


4 0 0 ? 9 . i 


5100 


DflTfl 


450? 9. 1 


5110 


DflTfl 


5 0 0 ? 9 . i 


5120 


DflTfl 


550? 9. 


5130 


DflTfl 


600 ? 9 . 1 


5140 


DflTfl 


6 5 O ? 9 . 


5150 


DflTfl 


700? 9. ‘ 


5160 


DflTfl 


750? 9. • 


5170 


DflTfl 


800 3 9 . ' 


5180 


DflTfl 


8 5 8 ? 9 . 


5190 


DflTfl 


900 ? 9 . ' 


5200 


DflTfl 


950 ? 9 . ‘ 


5210 


DflTfl 


1 000 3 9 


5220 


DflTfl 


1 050 ? 9 


5230 


DflTfl 


1 1 00 ? 9 


5240 


DflTfl 


1 150 3 9 


5250 


Dfl.Tfl 


1200? 9 


5260 


DflTfl 


1250 ? 9 


5270 


DflTfl 


1300? 9 


5280 


DflTfl 


1 350 ? 9 


5290 


DflTfl 


1 400 ? 9 


5300 


DflTfl 


1 450 ? 9 


5310 


DflTfl 


1500? 9 


5320 


DflTfl 


1550? 9 


5330 


DflTfl 


1 600 3 9 


5340 


DflTfl 


1 6 5 0 ? 9 


5350 


DflTfl 


1700? 9 


5 3 6 0 


DflTfl 


1 750 ? 9 


5370 


DflTfl 


1 800 ? 9 


5380 


DflTfl 


1 850 ? 9 


5390 


DflTfl 


1 900 ? 9 


5400 


DflTfl 


1950? 9 


5410 


DflTfl 


2080 ? 9 


5420 


END 





POWER 1 1A/CM+2" 
50? IE-50 

L 36E+ 1 3 ? 1 . 283E+13 
, 398E+13? 3. 083E+1: 
, 4 8 b E+13? 4 
, 7S2E+ 1 3 ? 5 
, 054E + 14.-5 
, 1 0 1 E + 1 4 ? 6 
. 1 4! 9 E + 1 4 ? b 
, 147E+14? 6 
, 157E+14? 6 



7 1 6 E + 1 y ? 



r' 1 6 E + 1 9 ? 

7 1 6 E + 1 9 ? 

7 1 6 E + 1 9 ? 



'16E+19? 

’16E+19? 



4 3 9 E + 1 3 
329E+13 
886E+13 
229E+13 
438E+13 
564E+13 
641 E+13 
687E+13 
715E+13 
3 4 E + 1 3 
742E+13 
748E+13 
751 E+13 
754E+13 
755E+13 
756E+13 
55 



’^E+l 



1 
1 
1 
1 

’ 16E + 19? 1 
1 
1 
1 
1 

7 1 6 E + 1 9 1 1 
7 1 6 E + 1 9 ? 1 
716E+19? 1 
716E+19? 1 
7 1 6E + 1 9 ? 1 
7 1 6 E + 1 9 ? 1 
716E+19? 1 
7 1 6 E + 1 y ? 1 
716E+19? 1 
7 1 6 E + 1 9 ? 1 
716E+19? 1 
, 7 1 6 E + 1 9 ? 1 



1 63 E + 14? 6 

1 6 7 E + 1 4 ? 6 
1 7 E + 1 4 ? b . 

171E+14- 6 
1 7 2 E + 1 4 ? 6 
1 7 3 E + 1 4 1 6 
1 7 3 E + 1 4 ? S 
1 7 3 1 + 1 4 j 6 

1 7 4 E + 1 4 j 6 
1 7 4 E + 1 4 ? 6 
. 175E+14? 6.757E+1 

. 175E+14? 6. 758E+1:: 
. 175E+14? 6. 758E+1:; 
. 1 7 6 E + 1 4 ? 6 . 7 5 9 E + 1 •. 
. 1 76E+ 1 4 ? 6 . 7 6E+ 1 3 
. 177E+14? 6. 76E+13 
. 1 7 8 E + 1 4 ? 6 . 6 7 1 E + 1 • 
. 1 7 8 E + 1 4 ? 6 . 7 6 2 E + 1 1 
. 1 79E+ 1 4 ? 6 . 7 63E + 1 •’ 
. 18E+14? 6. 765E+13 
. 1 8 1 E + 1 4 3 6 . 7 6 6 E + 1 • 
. 1 82E+ 1 4 3 6 . 7 68E + 1 • 
. 1 8 3 E + 1 4 ? 6 . 76 9 E + L 
,, 1 84E+ 1 4 ? 6 . 77 1 E+ 1 : 
. 1 8 6 E + 1 4 ? 6 « 7 7 3 E + 1 • 
. 1 8 7 E + 1 4 3 6 . 7 7 6 E + 1 • 
, 1 8 8 E + 1 4 ? 6 . 7 7 8 E + 1 • 
. 1 9 E + 1 4 ? 6 . 7 8 1 E + 1 3 
. 1 92E+ 1 4 ? 6 . 784E+ 1 : 
. 192E+14? 6. 787E+1: 
. 1 9 6 E + 1 4 ? 6 . 7 9 1 E + 1 ■ 
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TABLZ X. 



Output of the HP9330 for n-1 at Low Power. 



LOW POWER 10A/CM+2 



H 


- 


1 


IS 5 , 42 1 53E-05 


RT 


TIME 


0 


NSEC 


H 


- 


1 


IS 5 . 42 126E-05 


RT 


TIME 


50 


NSEC 


H 


- 


1 


IS 5.42097E-O5 


RT 


TIME 


100 


NSEC 


H 


- 


1 


IS 5.42O78E-05 


RT 


TIME 


150 


NSEC 


N 


- 


1 


IS 5.4206 7 E - 0 5 


RT 


TIME 


200 


NSEC 


N 


- 


1 


IS 5. 420S0E-05 


RT 


TIME 


250 


NSEC 


N 


- 


1 


IS 5 . 4205SE-05 


RT 


TIME 


300 


NSEC 


H 


- 


1 


IS 5 . 4205 3 E -05 


RT 


TIME 


350 


NSEC 


H 


- 


1 


IS 5 . 42052E-05 


RT 


TIME 


4 0 0 


NSEC 


N 


- 


1 


IS 5. 42O51E-05 


RT 


TIME 


450 


NSEC 


H 


- 


1 


IS 5 . 42O50E-05 


RT 


TIME 


500 


NSEC 


N 


- 


1 


IS 5.42 0 5 0 E - 0 5 


RT 


T I ME 


550 


NSEC 


N 


- 


1 


IS 5.42 0 50E- 0 5 


RT 


TIME 


600 


NSEC 


N 


- 


1 


IS 5.42 0 5 0 E - 0 5 


RT 


TIME 


650 


NSEC 


H 


- 


1 


IS 5 . 4205OE-05 


RT 


TIME 


7 O 0 


NSEC 


N 


- 


1 


IS 5.42049E-05 


RT 


TIME 


750 


NSEC 


N 


- 


1 


IS 5.42049E-05 


RT 


TIME 


800 


NSEC 


N 


- 


1 


IS 5.42049E-O5 


RT 


TIME 


850 


NSEC 


H 


- 


1 


IS 5 . 42049 E- 05 


RT 


TIME 


900 


NSEC 


H 


- 


1 


IS 5 . 42049E-05 


RT 


TIME 


950 


NSEC 


N 


- 


1 


IS 5 . 42049E-05 


RT 


TIME 


1 000 


NSEC 


H 


- 


1 


IS 5 . 42049E-05 


RT 


TIME 


1050 


NSEC 


N 


- 


1 


IS 5.42049E-05 


RT 


TIME 


1 100 


NSEC 


N 


- 


1 


IS 5 . 42049E-05 


RT 


TIME 


1150 


NSEC 


H 


- 


1 


IS 5 . 42049E-05 


RT 


TIME 


1 2 O 0 


NSEC 


N 


- 


1 


IS 5 . 42049E-85 


RT 


TIME 


1250 


NSEC 


H 


- 


1 


IS 5 . 42049E-O5 


RT 


TIME 


1 300 


NSEC 


N 


- 


1 


IS 5 . 42U49E- O 3 


RT 


TIME 


1350 


NSEC 


H 


- 


1 


IS 5.42049E-05 


RT 


TIME 


1400 


NSEC 


H 


- 


1 


IS 5.42049E-05 


RT 


TIME 


1450 


NSEC 


H 




1 


IS 5.42049E-05 


RT 


TIME 


1500 


NSEC 


N 


- 


1 


I S 5 . 42 0 4 9 E - 0 5 


RT 


TIME 


1550 


NSEC 


H 


- 


1 


IS 3.42 0 4 9 E — U 5 


RT 


TIME 


1 600 


NSEC 


H 


- 


1 


IS 5. 42049E-05 


RT 


TIME 


1650 


NSEC 


N 


- 


1 


Is 3.42 0 4 9 E “ O 5 


RT 


TIME 


1700 


NSEC 


H 


- 


1 


IS 5.42049E-05 


RT 


TIME 


1750 


NSEC 


H 


- 


1 


IS 5 . 42049E-05 


RT 


TIME 


1 800 


NSEC 


N 


- 


1 


IS 5 . 42O49E-05 


RT 


TIME 


1850 


NSEC 


H 


- 


1 


IS 5. 42048E-05 


RT 


TIME 


1900 


NSEC 


N 


- 


1 


IS 5 . 42048E-05 


RT 


TIME 


1 950 


NSEC 


N 


- 


1 


IS 5 . 42048E-05 


RT 


TIME 


2000 


NSEC 
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TABL2 XI. Output of the HP9330 for n-1 at High Power. 



HIGH POWER 1 lfl/CrtT-2 



N 




1 


IS 3. 42 153E-85 


Fi r 


TIME 


0 


MSEC 


N 


- 


1 


IS 5. 42123E-05 


AT 


TIME 


50 


NSEC 


H 


- 


1 


IS 5,42 0 9 1 E - 0 5 


AT 


TIME 


1 00 


MSEC 


N 


- 


1 


IS 5. 42071E-05 


AT 


TIME 


150 


MSEC 


H 


- 


1 


IS 5 . 42 9 58 E -05 


AT 


TIME 


200 


MSEC 


N 


- 


1 


IS 5.42051E-05 


AT 


TIME 


250 


MSEC 


H 




1 


IS 5 , 4 2 0 4 6 E “ 0 5 


AT 


TIME 


300 


MSEC 


H 


- 


1 


IS 5.42043E-05 


AT 


TIME 


350 


MSEC 


N 


- 


1 


IS 5.42042E-05 


AT 


TIME 


400 


MSEC- 


H 


- 


1 


IS 5 . 42041 E-05 


AT 


TIME 


450 


MSEC 


H 


- 


1 


IS 5.42040E-05 


AT 


TIME 


500 


MSEC 


N 


- 


1 


IS 5. 42040E-05 


AT 


TIME 


550 


MSEC 


H 


- 


1 


I S 5 . 4 2 0 3 9E- 0 5 


AT 


TIME 


600 


MSEC 


N 


- 


1 


IS 5 . 42039E-05 


AT 


TIME 


650 


MSEC 


H 


- 


1 


IS 5. 42039E-05 


AT 


TIME 


700 


MSEC 


H 


- 


1 


IS 5. 42039E-05 


AT 


TIME 


750 


MSEC 


N 


- 


1 


IS 5. 42039E-05 


AT 


TIME 


800 


MSEC 


N 


- 


1 


IS 5. 42039 E-05 


AT 


TIME 


850 


M S E C 


H 


- 


1 


IS 5 . 42039E-05 


AT 


TIME 


900 


MSEC 


N 


- 


1 


IS 5 . 42039E-05 


AT 


TIME 


950 


MSEC 


N 


- 


1 


IS 5. 42039E-05 


AT 


TIME 


1 000 


MSEC 


H 


- 


1 


IS 5 . 42039E-05 


AT 


TIME 


1 050 


MSEC 


H 


- 


1 


IS 5 . 42039E-05 


AT 


TIME 


1 100 


MSEC- 


N 


- 


1 


IS 5.42 0 3 9 E - 0 5 


AT 


TIME 


1150 


MSEC 


H 


■- 


1 


IS 5 . 42039E-05 


AT 


TIME 


1 200 


MSEC 


H 


- 


1 


IS 5.42 0 3 9 E-05 


AT 


TIME 


1250 


MSEC 


H 


- 


1 


I S j i 42 y 3 9 E — y 3 


AT 


TIME 


1 300 


MSEC 


N 


- 


1 


IS 5.42O39E-05 


AT 


TIME 


1350 


MSEC 


N 


- 


1 


IS 5. 42039 E-05 


AT 


TIME 


1400 


MSEC 


N 


- 


1 


IS 5.42 0 3 8 E - 0 5 


AT 


TIME 


1450 


NSEC 


H 


- 


1 


IS 5.42038E-05 


AT 


TIME 


1500 


MSEC 


N 


~ 


1 


IS 5. 42038E-05 


AT 


TIME 


1550 


MSEC 


H 


- 


1 


IS 5. 42038E-05 


AT 


TIME 


1 600 


MSEC 


N 


- 


1 


IS 5.42038E-05 


AT 


TIME 


1650 


MSEC 


N 


- 


1 


IS 5 . 42038E-05 


AT 


TIME 


1700 


MSEC 


H 


- 


1 


IS 5 . 42038E-05 


AT 


TIME 


1750 


MSEC 


H 


- 


1 


IS 5.42038E-85 


AT 


TIME 


1 8 0 0 


MSEC 


N 


- 


1 


IS 5. 42038 E-05 


AT 


TIME 


1850 


MSEC 


N 


- 


1 


IS 5 . 42037E-05 


AT 


TIME 


1900 


MSEC 


H 




1 


IS 5 . 42037E-05 


AT 


TIME 


1950 


MSEC 


H 




1 


IS 5 . 420 37 E-05 


AT 


TIME 


2000 


MSEC 
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VII. CALCULATION OF PHASE SHIFT 



3eam quality as previously mentioned depends on favor- 
able illumination of the exit aperture; favorable illumina- 
tion implies near constant phase. Therefore, calculation 
of the phase shift is important. To calculate the phase 
shift consider a laser cavity of length L as shown in Figure 



10 . 



CAVITY 



LASER OUTPUT 



L 



*1 



Figure 10. Geometry of Cavity 



Other symbols essential for the calculation are 



N - number of waves X - wavelength 

1 - denotes low power n - index of refraction 

2 - denotes high power v - frequency of the laser 



The number of waves in length L is given by 




( 59 ) 
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In equation (59) i = 1,2. Wavelength is related to index 
of refraction by 



‘ V* - nTv 



(60) 



Consider two regions in the cavity which are pumped at 
different power inputs. The difference in the number of 
waves in the two regions is 



II 



2 




(61) 



Substituting equation (60) into equation (61) yields 



No - N, = L (- 



n 0 v 



n^v 



) = 



Lv 



(n 2 -n 1 ) 



(62) 



Multiplying the right hand side of equation (62) above and 
below by n^ produces 

Lv n, n~-n, 

N 2 - n x - -5-i (-4^) (63) 



Multiplying equation (63) by 1/L and substituting equation 
(50) yields 



N 2 - N 1 



1 n 2 -n i 

h ( n i 



(64) 



The term on the left hand side is the phase shift per unit 
length of cavity. The phase shift at 1000 nanoseconds for 
a nonresonance wavelength is now calculated. 

A = 353 nm = 353.0E-7 cm (Figure 8) 
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From Tables X and XI the values of n are 



n 2 = 1 + 5.420392-05 



n x = 1 + 5.420492-05 



Substituting into equation (64) 




1 ,1 . 0000542039-1 . 0000542049 

353.0E-7 ^ 1.0000542049 



•) = -2 . 8327E-5 cm 



-1 



L 



The phase shift at 1000 nanoseconds for resonance is calcu- 
lated below. 

A = 351 . 62- 7 

The values of n are 

n 2 = 1 + (-3. 71283E-06) 
n 1 = 1 + (1 . 6602E-06) 

Substituting into equation (64) 

N 2“ N 1 _ 1 /. 999996287170-1. 00000 5602 N _ 0 ncAr , 0 „ 

L 351.6E-7 ' 1.000006602 ; cra 

The Strehl ratio, which is discussed by Born and Wolfe 
[21], provides a criterion for beam quality in terms of 
phase distortion. The ratio is 



I = I o exp [ - ( 2 tt 6 / A) 2 ] 



(65) 
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where I is beam intensity in the far field due to a phase 
distortion of 5/A. Based on equation (65), the criterion 
of 



6/A < 0.1 (66) 

is frequently invoked for satisfactory beam quality. Com- 
paring equations (65) and (61) one notes that 

N 2 - M x « 6/A (67) 



The maximum length of the laser cavity, L, for good 
beam quality is 



L = */A. 



N 2 “ N 1 



(n 2 -n 1 ) 

c-M) 



(63) 



Substitution of numerical values yields the following 
results : 

- maximum L for a laser operating at a nonresonant wave- 
length is 35 meters . 

- maximum L for a laser operating at a resonant wave- 
length is 0.1116 meters. 

For a nonresonant laser wavelength, a laser cavity with 
length, L, less than 35 meters will have 5/ A less than 0.1. 
However, if the laser wavelength is in resonance with neon, 
a length of laser cavity of only 0.1116 meters will cause 
6/ A = 0.1. R.ecall the power into the laser cavity was 
varied by 10 percent. 
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VIII. RESULTS 



For a power variation of 10 percent the phase shift at 
nonresonance is small = -0.0028 for a 100 cm 

cavity) and is dominated by the neon ground state whose 
population is approximately one million times larger than 
the other species. Hox;ever, in the resonance case, if the 
assumed values of R* are correct, the index of refraction 
can be greatl;/ magnified, resulting in a phase shift of 
-8.959E-3 cm ^ as compared to -2.3327E-5 cm ^ for nonreso- 
nance . 

It should be noted that the populations of the neon 
excited states were calculated using an arbitrary choice to 
determine the rate constants. However, the population of 
Ne is approximately six orders of magnitude larger than 
the neon excited states and thus dominates. As long as the 
laser is operated away from the 3s -* 4p resonant wavelengths, 
the phase shift will be negligibly small resulting in 
satisfactory beam quality. 
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APPENDIX A 

KINETIC REACTIONS AND RATE CONSTANTS FOR. XeF STMIIT.ATTOM 
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21 > UC19) X^* Xt : E r > XE2* ML 1.6CCF-32 SP5d 
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APPENDIX B 



CALCULATION OF RATE CONSTANT 

Consider a thin layer of active optical medium, with 

gain a, cm”\ in the laser cavity of length l, cm. The 

2 

change in intensity I, watts/cm , within the thin layer is 



AI - W - r in * a 1 1 



(69) 



The number of photons per unit time entering the layer of 
active optical medium is 



N. 



m 




( 70 ) 



A corresponding term for the number of photons leaving per 
unit time is 



I _ 
l, = _°'i£ 
out hv 



(71) 



The rate of production of photons within the thin layer is 

• • 

the difference between II . and N. which is 

out m 

| rate of pro- 
duction 
of photons 

Equation (72) neglects spontaneous emission. Combining 
equations (69) and (72) yields 



N -N. = A- (I -I.) 

out m hv out m 



(72) 
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production rate = — 



(73) 



From equation (19) the rate of production of photons by 
stimulated emission is 

production rate = k^ [hv] [KrF*]A£ (74) 



where the symbol [hv] means number of photons per unit volume. 
The gain coefficient also is equal to 

a = a [ KrF* ] (75) 



where a is the optical cross section for stimulated emission 
at the frequency v. The number density of photons [hv] is 
related to intensity by 



[hv] 



I 

chv 



(76) 



where c is the velocity of the photons; c, obviously, is the 
speed of an electromagnetic wave. 

Combining equations (73) to (76) gives 

h 5 = ca (77) 

3 

The dimensions of are cm /sec; the dimension of o is cm, 

2 

and c is cm /sec. From reference [2] the value of a is 
2 

5.021-16 cm , and c is 3.0E10 cm/sec. Substituting these 
values into equation (77) yields 



77 



= 1.5E-05 cm /sec 



This is the value for used in the computer program. 
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APPENDIX C 



CALCULATION OF PHOTON CONCENTRATION 
The concentration of photons is given by equation (76) 



as 


[hv] = I/chv 


Assuming , 


I = 10 koj/cm (saturation value) 
c = 3.0S10 cm/sec 
h = 6 . 6256E-34 J-sec 
A q = 243.0E-9 m 



the concentration of photons can be calculated. 



[hv] • 


10 .0E3 

(3.0E10) (6 . 6256E-34) (^off) 


[hv] > 


= 4 . 159E11/ cm 3 . 
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